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ABSTRACT 


The  problem  of  electromagnetic  radiation  and  scattering  from  perfectly 
conducting  bodies  of  revolution  of  arbitrary  shape  is  considered.  The  math¬ 
ematical  formulation  is  an  integro- differential  equation,  obtained  from  the 
potential  integrals  plus  boundary  conditions  at  the  body.  A  solution  is 
effected  by  the  metnod  of  moments,  and  the  results  are  expressed  in  terms 
of  generalized  network  parameters.  A  computer  program  for  computing  the 
generalized  impedance  matrix  of  an  arbitrary  body  of  revolution  is  included. 

The  expansion  functions  chosen  for  the  moment  solutions  are  harmonic 
in  0  (azimuth  angle)  and  subsectional  in  t  (contour  length  variable).  Because 
of  rotational  symmetry,  the  solution  becomes  a  Fourier  series  in  0,  each  term 
of  which  is  uncoupled  to  every  other  tern.  Hence,  the  problem  reduces  to  a 
set  of  independent  modes,  one  for  each  harmonic  term. 

Illustrative  computations  are  given  for  radiation  from  apertures  and 
plane-wave  scattering  from  bodies  of  revolution.  The  impedance  elements, 
currents,  radiation  patterns,  and  scattering  patterns  for  a  conducting  sphere 
are  computed  both  from  the  general  program  dnd  from  the  classical  eigenfunction 
solution.  The  agreement  obtained  serves  to  check  the  general  program.  Sim¬ 
ilar  computations  for  a  cone- sphere  illustrate  the  application  of  the  general 
program  to  problems  not  solvable  by  classical  methods. 
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I. _ INTRODUCTION 


Determination  of  the  behavior  of  a  conducting  body  in  a  known  impressed 

field  is  a  fundamental  problem  of  applied  electromagnetic  theory.  Specific 

examples  are  scattering  by  conducting  objects  and  radiation  from  conducting 

antennas.  A  general  procedure  for  treating  such  problems  is  given  by  the 

[  1  2] 

method  of  moments.  *  This  is  a  procedure  for  reducing  the  functional 
problem  to  a  finite  dimensional  matrix  problem.  In  electromagnetic  theory, 
the  resultant  matrices  can  be  identified  as  generalized  network  parameters . 
Once  the  generalized  impedance  matrix  of  a  body  is  known,  the  behavior  of 
that  body  for  arbitrary  excitation  is  easily  calculated  by  matrix  manipula¬ 
tions.  Furthermore,  3  effect  of  impedance  loading  can  be  accounted  for  by 
constraints  on  the  matrix  equations. 

The  problem  is  formulated  in  the  conventional  way  as  follows.  Let  g1 
denote  the  known  impressed  field  and  £S  the  scattered  field  due  to  currents 
on  the  body.  Then  the  total  field  £  is  the  sum  of  the  impressed  and 
scattered  fields,  that  is, 

e=Ei  +  £s  O', 


The  scattered  field  can  be  expressed  in  terms  of  a  vector  potential  &  and 
scalar  potential  <|>  as 


t  -  -jHfc  -  (2) 


where 


TnF  ds 

(3) 

0  T3T  ds 

(*0 

2 


Here  S  is  the  surface  of  the  conductor,  R  is  the  distance  from  a  source  point 
to  the  field  point,  j£  is  the  surface  current  on  S,  and  o  is  the  surface 
charge  on  S.  The  current  and  charge  are  related  by  the  equation  of  continuity 

v  '•  il  =  -jwo  (5) 

The  boundary  condition  requires  that  the  tangential  component  of  total  £ 
vanish  on  S.  Hence, 

-  <n  <6> 

where  the  subscript  tan  denotes  tangential  component  on  S.  The  problem  can 
now  be  stated  succinctly  as 

Ui)  =  (7) 

where  L  is  the  integro- differential  operator 

LU)  =  3t>]tan  (8) 

A  solution  of  (7)  gives  the  current  jJ  on  S.  Usually  we  are  interested  in 
some  functional  of  which  can  be  computed  once  £  is  known. 

This  report  considers  the  special  case  of  bodies  which  are  rotationally 
symmetric  about  an  axis.  Because  of  this  symmetry  a  Fourier  series  expansion 
ir.  the  angle  of  rotational  symmetry  reduces  the  problem  to  a  system  of  inde¬ 
pendent  modes  This  is  important  from  the  standpoint  of  computation,  because 
it  is  faster  to  invert  several  small  matrices  instead  of  one  large  one.  The 

particular  case  of  scattering  by  conducting  bodies  of  revolution  has  been 

[6] 

treated  previously  by  Andreasen.  His  solution  is  conceptually  similar  to 
ours,  but  differs  in  detail. 
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II.  METHOD  OF  SOLUTION 


To  effect  *»  '"'''"tier,  cf  (7)  u&t  one  mtuioa  01  moments,  which  is 

Cl  O  1 

closely  related  to  Galerkin's  method.  *  This  procedure  approximates  (7) 
by  a  matrix  equation,  which  can  then  be  inverted  by  known  algorithms.  The 
matrix  so  obtained  is  a  generalized  impedance  matrix  for  the  body.^  The 
excitation  of  the  body  is  represented  by  a  voltage  matrix,  and  the  resultant 
current  on  the  body  is  represented  by  a  current  matrix. 

For  the  method  of  moments,  let  the  inner  product  be  defined  as 


<&*!>  = 


K  •  Z  ds 


(9) 


where  if  and  *[  are  tangential  vectors  on  S.  A  set  of  expansion  functions 

(sL }  is  next  defined,  and  the  current  on  S  approximated  by 
J 

(10) 

j 


where  I,  are  constants  to  be  determined.  Equation  (10)  is  substituted  into 
0 

(7),  which,  because  of  the  linearity  of  L,  reduces  to 


Z^ 


can 


(ID 


A  set  of  testing  functions  (Jf.  }  is  defined,  and  the  inner  product  of  ( 11 ) 
with  each  is  taken,  llie  result  is 

Z  <*i’  v  ’  <si> ji> 

j 


(12) 


i  =  1,2,5,...  •  The  subscript  tan  has  been  dropped  from  g1  because  the  inner 
product  involves  only  tangential  components.  We  now  define  the  generalized 
network  matrices 


[Z]  =  [<  )£.,  >] 

(13) 

[Vj  =  [<  g1  >] 

(14) 

[I]  =  Cr] 

(15) 

and  rewrite  the  set  (12)  as 

[z]  [ I j  =  [v] 

(16) 

[z]  is  the  generalized  impedance  matrix,  and  [y]  =  [z]  J' 

admittance  matrix.  The  inverse  of  (l6) 

is  the  generalized 

[i]  =  [y3  tv] 

(17) 

gives  the  coefficients  I.  of  the  current  expansion  (ID ), 

J 

approximate  solution  to  the  problem. 

The  impedance  elements  of  (13)  are  explicitly 

and  hence  is  an 

zij  -  $  %  •  +  ds 

(18) 

where  we  have  used  (8)  and  (9)-  The  subscript  j  denotes  that  4-  and  $ . 

0  J 

are  the  potentials  due  to  4-  and  a..  If  the  two-dimensional  divergence 

o  J 

theorem  is  applied  to  the  vector  $W_.  on  a  closed  surface,  the  following 


identity  results: 
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4  ■  «t  fi(t)  e'Jn!* 

4,i  -  fjtt)  e-Jn0 


(21*) 


which  differ  from  (22)  only  in  the  sign  of  the  exponent.  The  are  orthogonal 
to  ^  m  /  n,  over  0  to  2n  on  <1,  and  also  to  (the  field  from  jJ^).  Hence, 
all  impedance  elements  are  zero  except  those  for  which  m  =  n,  and  each  mode 
can  be  treated  separately.  This  is  the  major  simplification  introduced  by 
the  rotational  symmetry  of  the  body. 

The  use  of  (22)  and  {2m)  to  evaluate  the  elements  of  (2*1  )  results  in 
the  partitioned  matrix  equation 


"rz«] 

Cz^ll 

C i l) 

"iv^ 

n 

n  $ 

t 

n 

n 

Lz4 

tz^j  ! 

[4 

n 

"  J 

n 

n 

cf  the  Z 

suDmatrices  are 

(Ztt).  = 
r.  'ij 

<i> 

> 

j 

(Z^).  .  = 

CW1., 

ij 

> 

~niJ 

•nj 

(Z***-). .  , 
n  'ij 

*T.i’ 

> 

j 

(Zf0)ii  = 

n  ij 

<K 9 

•*m’ 

> 

j 

(25) 


(26) 


the  elements  of  the  V  subrcatrices  are 
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<Vn>i  ’  <’  *> 


(27) 


and  the  elements  of  the  I  submatrices  are  the  coefficients  in  (23).  Note 
that,  for  N  terms  in  the  Fourier  series  on  0,  there  are  N  sets  of  matrix 
equations  (25). 

The  solution  to  (25)  can  also  be  written  in  partitioned  form  as 


r — 1 

< 

3  r+ 

UJ 

J 

-p 

1 _ 

1 

< — » 

[r] 

n 

[V5*] 
n  J 

ynt] 

[Y00] 
n  J 

r— 1 

ts.  a 

M 

_ !zLj 

(28) 


The  Y  submatrices  must  in  general  be  obtained  after  inversion  of  the  entire 
Z  matrix,  and  are  not  the  inverses  of  corresponding  Z  submatrices.  However, 
the  -n  mode  matrices  are  related  to  the  +n  mode  matrices  by 


LYtt] 

i7 

1 - J 

= 

[YU] 

[-Y^] 

-r. 

-n 

n 

n 

1 

* — * 

3  <r+ 

[Y00] 
-n  _ 

_ytj 

[Y0P  j 
n 

(29) 


Hence,  only  the  n  >  0  mode  matrices  need  be  inverted.  The  proof  of  (29) 
follows  from  the  fact  that  the  Z  matrices  satisfy  the  same  equation  (see 
next  section),  and  this  symmetry  survives  matrix  inversion. 

Finally,  for  an  explicit  solution  we  must  choose  the  f^(t),  the  t  expan¬ 
sion  functions.  It  is  known  that  subsectional  expansions,  using  pulse  functions 
or  triangle  functions,  give  we 11- conditioned  matrices. 


The  current  must  be 
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differentiated  to  obtain  the  charge;  hence  it  is  preferable  to  use  triangle 
functions.  Furthermore,  a  triangle  expansion  gives  a  piecewise- linear 

[l  2 

approximation  which  converges  about  twice  as  fast  as  a  step  approximation.  * 
Finally,  if  pjJ  is  expanded  in  triangles  instead  of  £,  the  divergence  formula 
and  the  treatment  of  the  end  points  of  t  become  simpler.  Hence,  we  chose 

f  (t)  =  i  T(t-t  )  (30) 

l  p  i 

where  T  is  the  triangle  function 

f  1  -  [t|,  |t|  <  1 

r(t)  =  i  (3i) 

1°  |t|  >  i 

When  using  these  functions,  distance  and  frequency  are  scaled  so  that  the  t^ 
are  one  unit  npert.  As  iu  Figure  1,  t  is  zero  at  the  lower  pole  and  N  at 
the  upper  pole.  There  are  (N-l)  triangle  functions  with  peaks  at  1,2,3- . -N-l. 

Numerical  evaluation  cf  the  impedance  elements  is  quite  difficult  be¬ 
cause  of  the  complexity  of  the  formulas.  An  approximate  evaluation  was 
obtained  by  approximating  each  triangle  function  by  four  steps.  The  potential 
integrals  due  to  each  step  are  evaluated  at  the  center  of  all  steps  by  a 
numerical  integration.  The  second  integration,  represented  by  (l8),  is  then 
a; uroxiraated  by  the  value  of  the  integrand  ot  the  center  of  each  step, 
summed  over  the  four  step  approximation  to  the  testing  triangle.  The  details 
of  the  solution  are  gi ren  in  the  next  section,  and  the  computer  program  ard 
instructions  for  using  it  are  given  in  Appendix  A. 
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u*  ’  =  sin  v  sin  v'  cos(0-0')  +  cos  v  cos  v' 

111  '  Uw  =  '"Sin  v*  sin(0-0' ) 

~  (37 

U0  '  Itj.  =  sin  v  sin  (0-0'  ) 

*  Up  =  cos (0-0’ ) 

Here  v  is  the  angle  between  the  t  direction  and  the  z  axis,  being  positive 
if  points  away  from  the  z  axis  and  negative  if  points  toward  it. 
Changing  (0-0' )  to  a  new  variable,  and  expressing  the  sine  and  cosine  terms 
of  (37)  as  exponentials,  one  0  integration  in  (32)  can  be  performed.  The 
remaining  0  integration  defines  the  Green's  function. ^ 


cos  n0 


where  Rq  is  given  by  (35)  with  0'  =  0.  With  f  given  by  (30),  the  resultant 
expressions  for  the  impedance  elements  (26)  are 
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(Z^).  . 
n  '13 


gn+l  f  gn-l 


(Zt0).  . 
n 


-M 

C  f  r  gn+l  gn-l 

=  \  dt'  l  dt  I  jHiT(t'-i)  T(t-j;(sin  v  sin  v*  ~ - ^ - - 

C  0 

+  COS  v  cos  v’  gj  +  T’(t'-i)  T’(t-j)gnJ 

rN  rN  g  -  g 

=  1  dt'  I  dt  £  -HiT(t'-i)  T(t-j)  sin  v*  n+— ^ ^ 


>in  v' 


VI  Bn-1 


(Z^).  . 
n  'ij 


iin  v - -- 


J  j 

0  0 


(Z^).  . 
n  'ij 


*  -JhPT'(t'-i>  8n] 

N  N 

(  dt'  ^  dt  j^+HiT(t’-i)  T ( t- j )  sin  v  — ^ ^ 

0  'J0 

i;  N 

^  dt'  j  dt  [jHiT(t’-i)  T(t-j)  e—— 


.»  gn+l  +  gn-l 


Jjjr  *(»•-!} 


Here  T'  is  the  derivative  of  a  triangle  function 


(  1,  -1  <  t  <  0 

T'(t)  =  \ -l,  0  <  t  <  1 

I 

L  M  >  1 


which  is  a  pulse  doublet.  Further  evaluitior  of  (.59)  is  done  by  numerical 
me  chods . 

The  integrations  of  (39)  involve  many  different  integrands,  and  to  re¬ 
duce  the  number  of  integrations  the  following  approximations  are  made.  For 
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tne  t  integration,  the  T  function  is  approximated  by  four  pulses  of  amplitudes 

1‘k>  ^h}  3/" k ’  l/k>  &s  shown  in  Figure  2(a).  The  derivative  of  T  is  repre¬ 
sented  exactly  by  four  pulses  of  amplitude  1,  1,  -1,  -i,  as  ohown  in 

Figure  2(b).  The  functions  p.  sin  v,  and  cos  v  are  assumed  constant  over 
each  pulse,  equal  to  tneir  values  at  the  midpoint  of  the  pulse.  For  the  f 
integration,  the  T  function  is  approximated  by  four  impulse  functions,  of 
strengths  1/8,  3/8,  3/8, 1/8,  as  shown  in  Figure  2(a).  The  derivative  of  T 
is  approximated  by  four  impulse  functions  of  strengths  1/2,  1/2,  -1/2,  -1/2, 
as  shown  in  Figure  2(b).  We  now  define  the  numbers 


CO 

H 

it 

.-•i 

t-3 
i — ’  — 

11 

ru 

a2  -  VS 

T’  =  1/2 

t3  =  5/8 

=  -1/2 

(M) 

co 

35* 

n 

-s- 

*3 

II 

» 

/o 

the  midpoints  of  the  pulses 

t  =  i  + 

+  _  ,  .  q  -  2.5 

0*2) 

P  2 

q  J  2 

and  the  pulse  Green's  functions 
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Then  jR^  -  R  |  <  1/4.  If  the  distance  between  peaks  of  the  triangle  functions 

is  a  tenth  of  a  wavelength,  then  kjE^  -  |  <  n/20,  and  exp  (-jkR^)  may  be 

represented  by  two  terms  of  its  Taylor  expansion  about  R  .  This  gives 

pq. 


M 


2  it 

u  =  "t; 

n  M 


Z-jkR 

cos  n0  e 
rm 

m=l' 

J  2 


I 


1  -  jk(R  -R  ) 

“ . ~5"  p--m- 


(«) 


For  the  final  integration,  we  approximate  t  by  a  straight  line  between 


=  0  +  ^~and  t  =  j  +  .  ?] 

of  (48)  can  be  simplified  by  completing  the  square.  'This  gives 


t  =  j  +  :i~and  t  =  j  +  .  Then  z  and  p  are  linear  in  t,  and  the  integrand 


M 


r 

Gn  ~  M 


t  +  1/4  ,—x — *- 

-jkR  °r  1"  jk(  yz+a  -  R  ) 

cos  n0m  e  1X1  \  dt  - 7^3- g-  ^ •  (49) 


m=l 


I. 


to  - 1/4 


+  d 


where 


t  =  |(z  -z  )  cos  v  +  (n  cos  0  -  0  )  sin  v  i 

o  P  q  9  VMp  rm  q  q1 


bo) 


2  2  2 
d  =  R  -  t 
pq  o 


(51) 


The  integration  of  (49)  is  now  straightforward,  and  one  obtain 


r-] 

<  si  n.--:  '  J 


M 


i  =  r,  /  cos  n0  f(0  ) 
n  M  /  rm  Mir 


m=l 


(52) 


where 


-jkR  f  t  +  VOd2 

L2(1+3V  l08rT^?'jkJ 


t,  =  t  -  T- 
1  0  4 


t_  =  t  +  y 
2  o  4 


(53) 


(54) 


This  completes  the  evaluation  of  the  impedances . 

Some  computational  details  will  now  be  discussed.  In  (54),  if  t  is 
larger  than  1/4  both  t^  and  t^  are  positive,  but  if  tQ  is  less  than  1/4,  t^ 
is  negative.  In  this  case  the  log  term  of  (53)  should  be  replaced  by 


(*2  + 


h?  +  a  )  (-tL  + 


42  +  d2) 


Gr  is  the  scalar  potential  due  to  a  charge  distribution  on  the  lateral  sur¬ 
face  of  a  frustum.  The  approximation  (46)  supposes  line  charges  at  0^.  If 
the  source  point  is  very  close  to  the  field  point,  the  granular  nature  of  the 
line  charges  may  be  too  restrictive.  When  the  source  and  field  points 

coincide  on  the  first  interval  0  <  0  <  -n,  the  line  charge  at  0  =  _£  is  spread 

-  r  -  W  2M 

out  into  a  uniform  surface  charge.  Furthermore,  if  the  surface  of  this  first 

Pq51 

interval  is  approximated  by  a  plane  rectangle  of  dimension  1/2  by  -g-  , 


(56) 


Now,  since 
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one  can  reduce  (56)  co 


A  more  elaborate  method  involving  a  localised  increase  in  the  density  of 
line  charges  was  attempted  but  wa?  abandoned  when  the  results  did  not  differ 
appreciably  from  (53). 
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IV.  COMPUTATION  OF  IMPEDANCES 


Equation  (45)  for  the  impedances  has  been  programmed  in  Fortran  IV  for 
execution  on  an  IBM  System/360  Model  50  computer.  Appendix  A  lists  the  com¬ 
puter  program  and  all  its  printed  output  for  the  first  mode  (n-1  in  (45)) 
for  a  sphere  of  radius  ^  and  raaius  to  wavelength  ratio  0.2.  Also  included 
in  Appendix  A  are  instructions  for  the  use  of  the  computer  program  along 
with  some  explanatory  notes  to  facilitate  possible  modifications.  In  the 
interest  of  simplicity,  the  program  is  written  tc  compute  one  mode  at  a 
time.  Input  data  consists  of  the  mode  number,  the  total  length  (an  integer 
determining  the  sizc  of  the  impedance  matrix  desired)  of  the  curve  that 
generates  the  surface  S,  ar.  integer  telling  how  many  times  to  subdivide  the 
0  axis  between  0  and  l8c^  for  the  Gn  of  (52),  the  wave  constant  and  some 
arrays  defining  the  angle  between  the  tangent  to  the  generating  curve  and 
the  z  axis,  the  cylindrical  coordinate  radius  p  of  the  generating  curve  and 
z  the  axial  distance  from  the  start  of  the  generating  curve  at  the  lower 
pole.  The  output  consists  of  some  of  the  of  (52),  the  impedance  matrix 
Z,  and  its  inverse  the  admittance  matrix  1 ,  as  explained  in  Appendix  A. 

To  check  the  accuracy  of  the  program,  the  impedance  for  current  on  a 
sphere  has  been  evaluated  by  spherical  mode  expansions.  The  details  of  the 
solution  are  given  in  Mautz’s  thesis.  The  final  result  for  any  mode 

with  e"  variation  on  a  sphere  of  radius  a  <  „ 
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where 


D  _  ^na‘~(n)(n+l)(n+m): 
mn  (2n+ 1 ) (n-m) ! 


(61) 


The  eigenvalues  are  given  in  terms  of  alternate  spherical  fessel  functions 
as 


=  tj  Jn(ka)  H^(ka) 


KnM  =  ^  Jn(ka)  Hi2)'  (ka) 


(62) 


The  b  and  c  coefficients  are 


nj 


r 


b  .  =  -2na  m  \  P“‘(cos  G)  f.(t)  dG 


:•  .  =  2~a‘-  f  - 

nJ  J 


11  sr-coc  G  " 

n 


(6;>) 


S3 n.  G  f  (t )  dw 


t 

where  f.(t)  is  tne  j'*'  expansion  function  as  defines  in  gen-rn  1  ;v  ;2  ;  and 
J 

m  particular  by  The  F^‘  in  (6*)  are  the  associated  Legendre  p^Voir.ials 

The  evaluation  of  '  o )  was  uor.e  nur-e  ricaiiv.. 
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Table  1  shows  representative  computations  for  a  sphere  of  radius  a  =  0.2>., 
using  9  equispaced  triangles.  In  other  words,  the  expansion  functions  are 
(90)  with  the  t^  =  1,2, ...,9,  and  the  z  axis  cuts  the  sphere  at  t,  =  3  anu  10. 
The  M  appearing  in  (52)  is  twenty.-  The  mode  chosen  is  n  =  1  in  (h5),  and 
correspondingly,  m  =  1  in  (60).  Column  1  of  Table  1  gives  the  computation 
using  the  general  solution  (U5)  of  Section  III.  Column  2  gives  the  solution 
using  IS  terms  of  the  infinite  series  (60).  Column  3  uses  9  terms  of  the 
series  (60).  Notice  that  the  real  parts  of  and  Z^  and  the  imaginary 
part  of  Z  are  identical  for  the  l8  term  and  9  term  solution  of  (60) 
and  are  only  slightly  different  from  those  of  (*+5)-  The  imaginary  parts  of 
Z4t  and  Z^  and  the  real  parts  of  and  Z^  of  (60)  are  converging  very 
slowly  in  the  modal  solution,  but  are  converging  towards  the  corresponding 
parts  of  the  general  solution.  Z^  is  supposed  to  be  zero,  and  the  non  zero 
values  of  Z^  in  Table  1  are  so  small  that  all  the  significant  figures  are 
probably  meaningless.  On  the  basis  of  comparisons  such  as  Table  1,  it 
appears  that  our  general  program  is  sufficiently  accurate  for  most  purposes. 


2.5 


V.  MEASUREMENT  MATRICES 


Any  linear  measurement  of  the  field  from  a  current  j[  on  a  body  S  can 
be  expressed  as  a  linear  functional  of  that  is 


measurement  = 


(64) 


where  JSr  is  a  known  function.  If  the  excitation-measurement 
as  a  two- port  system,  it  can  be  shown  that  £  in  (61)  is  the 

t  3  1 

the  measurement  port  is  excited.  For  a  moment  solution, 

given  by  a  superposition  ,£  =  2  I.,£.  ,  and  (61)  reduces  to 

3  3 


system  is  viewed 
field  on  3  when 
the  current  is 


measurement  =  [r]  [ij 


(65) 


where  [i]  is  the  matrix  (15)  and  [R]  is  a  measurement  row  matrix 

[R]  =  [<jL,  Er>]  .(66) 

J 

Note  the  similarity  of  [R]  to  the  excitation  matrix  defined  by  (ll).  If 
the  matrix  solution  (IT)  is  now  substituted  into  (65),  one  has 

measurement  =  [r][y][v]  {67) 

The  symmetry  of  (67)  with  respect  to  [R]  and  [v]  reflects  the  reciprocity 
theorem  of  electromagnetic  theory. 

For  bodies  of  revolution,  the  expansion  for  £  can  be  separated  into  t 
and  0  directed  components,  according  to  (23).  It  is  then  convenient  to 
partition  [r]  into  t  and  0  component  terms  as 


2k 


<Sn>i  ‘ 

(66) 

‘  <4-  ®r> 


The  analogous  partition  for  excitation  [V]  is  given  by  (?7)-  Now  one  can  re¬ 
write  (67)  in  partitioned  form  as 


measurement  =  £[r^]  [R^]  J 

J ' 

7v‘!_ 

1 — l 

"Sk 

L 

,tVn7 

(69) 


wnere  the  Y  submatrices  are  obtained  after  the  Z  matrix  is  inverted,  and  are 

not  the  inverses  of  the  corresponding  Z  submatrices. 

An  important  special  case  is  that  of  radiation  field  measurement.  It 

[  8  ] 

has  been  shown  that  the  radiation  field  from  currents  i  on  S  is  given  by 

£  '  31  =  e'jkr  [R]  [I]  (70) 


where  the  elements  of  [R]  are  given  by  (66)  with 


Er  =  u  e‘&  '  ^ 


(71) 


This  is  a  unit  plane  wave  with  polarization  vector  u  and  propagation  vector  k.. 
An  arbitrary  plane  wave  is  a  superposition  of  two  orthogonal  components,  say 
Eg  and  E^.  Hence,  one  can  treat  the  general  case  as  two  applications  of  (71)> 
one  for  ji  =  and  the  other  for  y.  =  u^.-  To  distinguish  between  the  two  cases 


let 
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(77) 
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where 

J  =  J  (kp  sin  8  ) 
n  n'  H  r 

Similarly,  to  evaluate  (73),  let 

jk(p  sin  8  cos  0  +  z  cos  8  ) 

n*  r  2?  I* 

The  required  dot-products  are  then 

*  ^  =  sin  v  sin  0 

^  •  u^  =  cos  0 

and  the  0  integrations  in  (73)  can  be  performed.  The  result  is 


(78) 


(79) 


(80) 


(B*  ),  -  2»jn+l 

n  1 


I 


N 


jkz  cos  8  J  .,+J  .. 

^  /  j_  \  _  r  .  n+1  n-1 

dt  p  fi(t)  e  sin  v - — - 


0 

N 


(Bl) 


).  -  2,Jn+1 


I 


dt  p  fjU)  e 


jkz  cos  8  J  ,,-J 

r  n+1  n-1 


where  (78)  again  applies.  For  computation,  the  p  f  (t)  in  ( 7 rJ)  and  (8l) 
were  the  triangle  functions  (31)-  The  remaining  integrals  must  be  evaluated 
numerically  for  bodies  of  arbitrary  contour.  In  general,  R^6  and  are 
even  in  n,  while  and  R^  are  odd  in  n.  As  shown  by  (22)  and  (2J+),  the 
excitation  matrix  [v]  differs  from  the  measurement  matrix  [R]  only  by  the  sign 
of  n.  Hence,  for  plane-wave  excitation  of  the  body, 


(VM).  =(RPq). 
n  '1  v  -n'l 


(82) 
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where  pq  represents  t9,  09,  t0,  or  00.  Equation  (82)  means  that  the  V  are 

given  by  (77)  and  (8l)  with  n  replaced  by  -n  and  9  by  9  ,  and  does  not  imply 

r  u 

equality  of  excitation  and  measurement. 

For  axial  incidence  of  plane  waves,  equations  (77)  and  (8l)  reduce 
considerably.  In  this  case  9^  =  n  and  all  Rn  =  0  except 

r  rN 

(R^9)i  =  (R^p^  =  -n  j  lit  p  f.(t)  e  ^kz  sin  v 

0 

N 

(R^0)i  =  -(R^p)p  =  j*  J  dt  p  f,(t)  e“^kz 

0 

(R^).  =  -(R^).  =  -j(Rf  l  (83) 

(Rf ).  .  (S?{).  =  -j(Rf ). 

Hence,  only  the  n  =  -1  and  n  =  +1  modes  are  excited  by  axially  incident  plane 


waves . 
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VI .  APERTURE  ANTEMAS 


An  aperture  antenna  is  a  radiating  system  consisting  of  a  conducting 

body  with  apertures  through  which  electromagnetic  energy  is  supplied.  For 

analysis  purposes,  it  is  assumed  that  the  tangential  electric  field  is 

known  over  the  aperture,  It  is,  of  course,  zero  over  the  conducting  body. 

s  i 

This  known  field  g^_  corresponds  to  E^.  =  -g^  in  the  general  analysis  of 
Section  I.  Figure  3  represents  the  general  problem,  showing  the  body  S, 
the  aperture,  and  radius  vector  to  a  distant  field  point. 

In  terms  of  a  moment  solution,  the  excitation  is  the  known  g^  in  the 
aperture,  from  which  the  excitation  matrix  may  be  computed  by  (l4)  by  setting 
if  =  -g^,.  In  general,  this  results  in  an  excitation  of  many  e^n^  modes,  each 
of  which  can  be  treated  separately.  The  current  on  S  for  each  mode  is  found 
by  inverting  (25),  and  the  total  current  is  the  sum  of  the  modal  currents, 
that  is 

n 

where  n denotes  the  mode.  Similarly,  the  radiation  field  for  each  mode  is 
found  from  (70),  and  the  total  radiation  field  is  the  sum  of  the  modal 
fields,  that  is 

E  =  ^  Jkr  [R  ]  [y  ]  [V  ]  (85) 

u  4nr  _  n  n  n 

n 

where  u  =  ©  or  0,  and  n  again  denotes  the  mode. 

For  simplicity,  illustrative  computations  have  been  made  only  for  the 
case  of  excitation  independent  of  0,  which  excites  only  the  n  =  0  mode. 


in  aperture 


I 

I 


Figure  3*  An  aperture  in  a  conducting  body  revolution 
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This  results  in  a  furtr.er  simplification  in  that  [h^J  =  [z^]  =  0,  that  is, 
there  ^s  no  coupling  between  t»-e  t-uireeted  currents  and  the  ^-directed 
currents*  From  symmetry  considerations,  an  excitation  having  only  an  E 


component  produces  only  a  J  current,  and  an  excitation  having  only  an  E 

0  ( 

component  produces  only  a  current.  Hence,  +he  problem  reduces  to  two 
independent  ones, 


’0 


J.  »  Ejiltvn  rv:j 

t  0  0  j 


j0  -  Uj?J  [#1  uvgi 


(66) 


for  the  cut  rent,  and 


E 

b 


„  iJiik  e"-Xr 


]  [YU3  Tv1] 
J  u0  0 


(37) 


for  the  radiation  field.  The  lY^]  and  [Y^/1  in  (86)  and  (87)  are  the  inverses 

V  V 

of  [Z^j  and  [zfr],  respectively.; 

To  check  the  program,  radiation  from  an  equatoria Liy  slotted  conducting 
sphere  was  computed  both  by  the  classical  eigenfunction  expansion  and  by  the 
general  program.  Figure  t  shows  the  current  density  on  a  sphere  of  radius 
a  =  0.2a  excited  by  an  E_  =  6(©-jt/2)/a  in  the  equatorial  slot.  This  gives 
rise  to  V  “7  .  ~  2jt  if  t..  is  on  the  equator,  and  rero  otherwise.;  The  eigenfunction 

'  w  O 

solution  was  computed  using  15  terms  of  the  Fourier  series.  The  real  part  of 
the  current  converged  very  rapidly,  and  can  be  considered  as  exact.  The 
imaginary  part  of  the  current  does  not  converge  at  the  soi-iree,  and  converges 


a  conducting  sphere  of  radius  0.'2\,  excited  by 
equatorial  slot. 


Figure  5-  Current  density  on  a  conducting  sphere  of  r-t.lius 
impulsive  in  an  equatorial  slot. 
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slowly  elsewhere,  and  hence  the  imaginary  part  computed  from  13  terms  of  the 
Fourier  series  is  not  exact.  As  can  be  seen  in  Figure  1,  the  N  =  10  solution 
using  the  general  program  is  ir.-  good  agreement  with  the  eigenfunction  solution. 
Computation  of  radiation  field  patterns  for  the  two  solutions  agreed  to  within 
less  chan  one  percent  accuracy,  and  hence  are  not  shown. 

Figure  5  shows  the  current  density  on  a  sphere  of  radius  a  =  0.2A  excited 

by  an  E^  =  -6(9-n/2)/a  in  the  equatorial  slot.  This  gives  rise  to  =  2n 

when  t,  is  on  the  equation,  and  zero  otherwise.  Again  the  eigenfunction 
J 

solution  was  computed  using  rS  terms  of  the  Fourier  series,  and  the  real  part 
of  the  current  converged  very  rapidly.  However,  the  imaginary  part  of  the 
current  does  not  converge.  The  real  part  of  the  current  computed  from  the 
N  -  10  case  of  the  general  program  again  agrees  well  with  the  eigenfunction 
solution.  It  is  difficult  t,o  say  anything  about  the  imaginary  part  because 
of  the  divergent  nature  of  the  solution.  Computation  of  the  radiation  field 
patterns  for  the  tve  solutions  again  agreed  to  within  a  few  j  vreent. 

Evidently  the  radiation  pattern  docs  not  depend  strongly  on  the  non- converging 
imaginary  part  of  the  current. 

An  an  example  of  a  body  of  more  complicated  shape,  computations  were 
made  for  apertur*'  excitation  of  a  cone-sphere.  This  body  is  formed  by  wedging 
a  conducting  sphere  into  a  conducting  cone,  resulting  in  a  shape  similar  to 
that  of  an  ice  cream  cone.  The  radius  cf  the  sphere  is  a  =  Z.2\,  and  the 


cone  half  angle  is  ItT.  The  aperture  is  e.  narrow  slot  near  the  cone-to-sphere 
junction.  Figure  6  shows  the  current  density  wh«.  the  excitation  is  an 
impulsive  E  In  th^  slot.  An  N  =  20  solution  is  shown,  but  smaller  values  of 

v 

N  give  reasonably  accurate  solutions. 


5^ 


Figure  6.  Current  density  on  a  conducting  cone- sphere  excited  by  one  volt 
across  a  slot  at  the  cone- to- sphere  junction. 
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Current  density  on  a  conducting  cone-sphere  excited  by  an  impulcive 
in  a  slot  at  the  cone-to-sphere  .junction.. 


Figure  7  shews  the  same  cone-sphere  excited  by  an  impulsive  in  the 
slot.  Again  an  N  =  20  solution  is  shown,  but  usable  results  could  be  obtained 
with  a  smaller  N.  Note  the  extreme  osci  .latory  nature  of  the  imaginary  part 
of  the  current.  Theoretically,  the  imaginary  part  does  not  converge,  as  dis¬ 
cussed  previously  for  the  sphere. 

Figures  8  and  9  show  the  power  gain  patterns  for  the  two  cone- sphere 
excitations  discussed  above.  For  a  given  polarization,  the  power  gain  normal¬ 
ized  to  omnidirectional  radiation  is'--^ 


4*r2|Ej2 


T)  F 


where  tj  =  -/u/e  and  is  the  power  radiated.  From  (70)  and  (17),  the 

numerator  of  (88)  can  be  expressed  as 


4*r2  |EJ2  =  jL  ^  [R]  [y]  [V]  2 


In  general,  the  power  radiated  is  given  by 


P  =  -Re 
rad 


g*  •  I  ds  =  Re  [V]  [Y]  [V] 


where  the  elements  of  the  row  matrix  [ V 3  are 


S 


For  the  excitations  used  in  Figures  8  and  9>  only  one  is  nonzero,  and 

A 

=  2jt.  Hence  (90)  reduces  to  a  single  element  of  the  [Y]  matrix.  In 


Figure  8.  Power  gain  pattern  for  radiation  from  a  conducting  cone-«phere  excited  by  one  volt  across 
slot  at  the  cone- to* sphere  junction. 


particular,  for  Figure  8, 
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\  is  the  point  of  excitation.  Similarly,  for  Figure  9, 
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again  t^  is  the  point  of  excitation. 


VII.  PLANE-WAVE  SCATTERING 
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=  bn  j  [ar3  i s j  [a'j  1 


(102) 


where  [s]  is  given  by  (9b).  Equation  (9?)  is  a  special  'ase  of  (102). 

For  simplicity,  only  excitation  by  axially  incident  plane  waves  it 
considered  for  computation.  Let 


E1  =  -u  e-^2 
~  ~x 


(103) 


which  can  be  considered  a  ©-polarized  plane  wave  from  a  transmitter  at 

coordinates  9^  =  n,  =  0.  The  excitation  matrix  elements  are  thus  given 

by  the  first  two  equations  of  (83)  witn  R  replaced  by  V,  n  by  -n,  and  9  by  9  . 

r  x 

Only  the  n  =  +1  modes  are  present  in  the  excitation,  giving  rise  to  only 
the  same  two  modes  in  the  current.  The  symmetry  relationships  (29)  and  (32) 
can  now  be  used  to  obtain  the  simpler  form 


+■  $ 
si  =  cos  0  -r  i sin  0 


(104) 


where 


N-l 

jt  ■ 2 


dpi  f,(p) 


(105) 


*«2_ 


and  f.(p)  =  T(t-t. )/p.  Hence,  the  current  is  completely  specified  by  itr 
n  =  1  mode  coefficients.; 

To  check  the  program,  the  current  on  a  sphere  excited  by  the  plane  wave 
(1C3)  was  computed  both  by  the  classical  eigenfunction  method  and  by  the 
general  computer  program.  Figure  11  shows  the  magnitudes  of  J°  and  of  (10b) 


•nuu.-t  '.n, a ;  ii-  ]•>  of  rsiiiuG  C  . tx-ited  oy 


for  n  sphere  'f  rr-alus  a  =  G-2-%  using  N  -  10  for  the  general  program.  The  • 
t  ig-nfunctior.  cvqanslcn  converges  very  rapidly  and  can  be  considered  as  exact. 
The  car rent  density  is  normalised  with  respect  to  the  incident  magnetic  fi^ld 
intensity,  giving  a  dimensionless  quantity  which  is  2  for  a  plane  conductor 
covering  the  ?.  ~  C  plane,  The  radar  scattering  patterns  computed  by  the 
eigenfunction  method  and  by  the  N  =  10  general  solution  were  almost  identical, 
and  are  not  shown. 

Figure  12  snows  the  current  on  a  cone-sphere  excited  by  a  plane  wave 
axially  incident  on  the  tip..  The  radius  of  the  sphere  is  again  C.2A,  which 
is  the  same  as  used  in  the  antenna  examples.  Again  j''  and  Y  of  (iCw), 
normalized  with  respect  to  K“,  are  plotted.  To  illustrate  convergence,  both 
the  N  =  20  (circles)  and  the  N  =  30  (solid  line)  solutions  are  plotted.  The 

*■  H. 

J“  current  converges  fairly  rapidly,  but  the  -Y  current  converges  more  slowly. 

$ 

Even  tne  N  -  30  solution  cannot  be  considered  accurate  for  J  .  Apparently 
the  rate  of  convergence  is  affected  by  the  sharp  tip  of  the  cone.  For  better 
convergence  one  could  include  a  term  in  the  current  expansion  to  properly 
represent  the  singularity  at  the  tip.  However,  radar  scattering  patterns 
are  .insensitive  to  small  oscillations  in  the  currents,  and  more  accurate 
solutions  are  not  needed  in  most  cases. 

Figure  13  shows  the  current  on  the  same  cone-sphere  excited  by  a  plane- 
wave  axially  incident  on  the  spherical  end,  that  is,  the  end  opposite  to  the 
point  The  above  comments  concerning  convergence  apply  also  to  this  excitation. 

Figures  ll  and  15  show  radar  cross  section  patterns  for  tne  cone-sphere 
excitations  corresponding  to  Figurs  12  and  13,  respectively .  The  riots  were 
made  using  the  N  =  2C  solution,  but  the  corresponding  N  =  solution  gave 
almost  identical  results.  The  only  difference  occurred  ir.  the  vicinity  of 
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Current  density  on  a  conducting  cone-sphere  (cone  angle  20 ,  sphere  radius 
0.2 \)  excited  by  a  uniform  plane  wave  axially  incident  or.  the  sphere. 


forward  scattering  where  a  change  in  a  of  the  order  of  a  few  percent  was 

obtained.  This  region  is  most  sensitive  to  inaccuracies  because  all  parts 

tend  to  radiate  in  phase  in  the  forward  direction.  All  plots  are  normalized 

to  k  to  make  them  dependent  only  on  length/ wave length  dimensions,  not  on 

0©  2 

absolute  dimensions.  The  curves  labeled  a  /k  are  in  the  0  =  0  plane,  and 
are  the  radar  cross  sections  measured  by  a  0-polarized  receiver.  The  curves 
labeled  0^®  are  in  the  0  =  jt/2  plane,  and  are  the  radar  cross  sections 
measured  by  a  0-polarized  receiver. 

Because  only  the  n  =  1  and  n  =  -1  modes  are  present  in  the  excitation, 
these  are  the  only  modes  present  in  the  scattered  field.  From  (95)  and.  the 
symmetry  relationships  (29)  and  (32)  one  can  show  that 


where  and  s^S  are  the  n  =  1  modal  solutions  evaluated  at  0r  =  0.  Hence, 
from  (99)  and  (106)  it  follows  that  the  radar  cross  sections  are  given  by 


(107) 


Hence,  the  graphs  of  Figures  14  and  15  are  proportional  to  the  coefficients 


lettering  cross  section  for  conducting  cone-i 
ed  by  a  uniform  pLane  wave  axially  incident 


Figure  1 5.  Bistatic  scattering  cross  section  for  a  conducting  cone-sphere  (cone  angle  20°,  sphere  radius 
0.2\)  excited  by  a  uniform  plane  wave  axially  incident  on  the  sphere,  a®®  is  in  the  plane 


hi 
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VIII.  DISCUSSION 

This  report  develops  a  solution  for  the  electromagnetic  behavior  of  con¬ 
ducting  bodies  of  revolution  with  arbitrary  excitation  in  terms  of  generalized 

[  5  ] 

network  parameters.  basically,  ^his  involves  the  application  of  a  moment 

method  to  the  superposition  integral  representation  of  the  problem.  Because 
of  the  rotational  symmetry,  a  Fourier  series  expansion  on  0  is  used.  To 
facilitate  the  handling  of  arbitrary  body  contours,  a  subsectional  expansion 
on  t  is  used.  The  simplest  subsectional  function  is  the  pulse  function,  but 
it  is  not  in  the  domain  of  the  original  operator  because  derivatives  must  be 
taken.  The  derivative  can  be  replaced  by  a  finite  difference  approximation, 
thereby  approximating  the  operator.  However,  it  was  found  that  difference 
approximations  did  not  yield  as  good  accuracy  as  desired.  The  next  simplest 
subsectional  function  is  the  triangle  function  ,  and  this  was  used  for  the 
general  program.  Its  derivative  is  a  doublet  pulse,  and  hence  no  approxima¬ 
tion  of  the  operator  is  necessary. 

Three  types  of  singularities  are  encountered  in  the  general  solution: 

(a)  a  singularity  in  the  coordinate  system  along  the  z  axis,  (b)  singularities 
of  curvature  in  the  body  contour,  such  as  at  a  cone  tip,  and  (c)  singularities 
in  the  excitation.  The  first  type  introduces  difficulties  in  the  numerical 
evaluation  of  the  impedance  elements,  but  can  be  overcome  by  careful  analysis. 
The  second  type  could  be  taken  care  of  by  special  subroutines  at  points  of 
discontinuity  in  body  curvature,  but  would  be  difficult  to  implement.  No 
attempt  has  been  made  to  do  this  in  the  general  program.  The  third  type  of 
singularity  is  impractical  to  incluie  in  a  general  solution,  because  then  the 


impedance  matrix  would  depend  on  the  excitation  instead  of  being  a  function 
of  the  body  geometry  alone.  It  could,  however,  be  incorporated  into  particular 
solutions.  For  example,  the  singularity  i.  the  current  at  feed  point  of  a 
slot-fed  antenna  can  be  accounted  for  by  including  the  proper  singularity  in 
the  expansion  functions. 

If  only  the  radiation  and  scattering  patterns  are  wanted,  but  not  the 
current,  then  less  accuracy  is  required  in  the  solution.  This  is  because 
such  patterns  are  continuous  linear  functionals  of  the  current,  and,  as  shown 
by  the  calculus  of  variations,  are  accurate  to  the  second  order  when  the 
currents  are  accurate  to  the  first  order..  For  example,  if  the  current  has  an 
error  of  the  order  of  ten  percent,  the  radiation  field  has  an  error  of  the 
order  of  only  one  percent.  Furthermore,  the  radiation  depends  primarily  on 
che  current  moment  pj  instead  of  the  current  itself.  Hence,  inaccuracies  at 
the  poles  (p  =  0)  of  the  body  have  little  effect  on  the  radiation  pattern. 
Finally,  rapid  oscillations  in  the  current  about  the  correct  value  also  have 
little  effect  on  the  radiation  field.  They  contribute  mainly  to  near-field 
stored  energy. 

The  generalized  network  parameters  of  a  body  are  basically  a  matrix 
approximation  to  the  functional  operator  equation  for  the  problem.  Hence,  so 
long  as  the  computational  approximations  are  valid,  one  can  compute  the 
response  of  the  body  to  arbitrary  electromagnetic  excitation.  Aperture 
antennas  (Section  VI)  and  plane- wave  scattering  (Section  VII )  are  but  two 
examples.  Another  problem  of  interest  is  excitation  by  current  dipoles  in 
the  vicinity  of  the  body.  Still  more  generally,  both  the  excitation  and 
measurement  can  be  taken  in  the  near- zone  if  desired.  Radiation  by  loaded 
aperture  antennas  and  scattering  by  loaded  scattering  antennas  can  also  be 
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treated  by  known  methods .  , 

Other  problems  can  be  treated  by  modifications  or  extensions  of  the 

theory.  For  example,  the  internal  resonances  of  a  cavity  of  revolution  can 

be  found  by  a  procedure  similar  to  that  used  for  waveguides  of  arbitrary  cross 

[12] 

section.  The  problem  of  scattering  by  a  homogeneous  dielectric  or  magnetic 

body  of  revolution  can  be  treated  in  terms  of  equivalent  electric  and  magnetic 
surface  currents.  Hence,  an  extension  to  include  both  electric  and  magnetic 
currents  could  be  used  to  solve  such  problems. 


APPENDIX  A 


THE  COMPUTER  PROGRAM 


This  appendix  provides  a  description,  in  terms  of  Fortran  IV  language, 
of  the  computer  program  that  finds  the  generalized  impedance  and  admittance 
matrices  for  a  conducting  body  of  revolution  with  surface  S.  Section  A-l 
gives  instructions  for  using  the  computer  program  while  Section  A-2  gives  some 
explanatory  notes  to  facilitate  possible  modifications.  Section  A- 3  lists  the 
computer  program  along  with  the  data  for  the  first  mode  (n=l  in  (45))  of  a 
sphere  of  radius  10/ir  and  radius  to  wavelength  ratio  0.2.  Section  A-4  gives 
the  printed  output  for  this  same  data. 

A-l  Instructions  for  using  the  computer  program 

The  program,  written  in  Fortran  IV  level  G  and  tested  on  an  IBM  System/360 
model  50  computer,  accepts  punched  card  input  data  and  outputs  the  general!  zer. 
impedance  and  admittance  matrices  via  both  the  printed  line  and  a  direct 
access  storage  device. 

The  direct  access  storage  device  was  assigned  the  data  set  reference 
number  4.  Statement  number  1  early  in  the  main  program 


1  REWIND  4 


returns  the  data  set  number  4  to  the  first  record  if  it  is  not  already  there. 
Each  unformated  WRITE  (4)  statement  defines  a  new  record  on  data  set  number  4. 
If  the  first  n  records  have  been  written  previously,  it  may  De  desirable  to 
skip  them  and  start  writing  on  the  (n+1)  record.  This  is  effected  by  following 
statement  1  with  n  identical  READ  (4)  statements.  There  are  no  such  READ  (4) 
statements  in  the  listing  of  the  program  in  Section  A- 5. 
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All  the  punched  card  data  are  read  early  in  the  main  program  according 
to  the  statements 

50  READ  (1,51,END=52)  NN,  N,  NIKI,  BK 

51  FORMAT  '313, Elk. 7) 

KG  =  2*N 

61  READ  (1,53)(VS(J),  J  =  1,KG) 

62  READ  (1,53){ZS(J),  J  =  1,  KG) 

63  READ  (1,53)(R(J),  J  =  1,  KG) 

53  FORMAT  (10F8.4) 


in  that  order.  According  to  statement  50,  several  sets  of  data  may  be  pro¬ 
cessed  with  control  going  to  statement  52  when  the  data  are  exhausted. 
Statement  52  is  a  stop. 

MW  =  the  mode  number  the  same  as  n  appearing  in  (45 ) .  MN  must  be 

either  zero  or  a  positive  integer.  If  NN  is  large,  NIKI  must  also 
be  large  for  the  numerical  integration  (46)  to  be  accurate. 

N  =  total  length  of  the  generating  curve.  N  is  a  positive  integer. 
There  will  be  N-l  expansion  functions  for  the  t  directed  current 
and  N-l  expansion  functions  for  the  0  directed  current.  Conse¬ 
quently,  the  Z  end  Y  matrices  will  be  2*(N-l)  by  2*(N-l).  For  N 
larger  than  20,  present  dimension  statements  must  be  altered. 

NIKI  =  the  number  of  equal  subdivisions  of  the  0  axis  from  0°  to  l80°. 
NIKI  corresponds  to  the  positive  integer  M  appearing  in  (46). 


NIKI  must  be  large  enough  so  that 


NPHI 


*  NN  is  less  than  a  radian. 


The  maximum  electrical  length  corresponding  to  an  excursion 


NHil 


in 
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0  should  also  be  less  than  a  radian.  If  NPHI  is  larger  than  40, 
present  dimension  statements  must  be  altered. 
dK  =  propagation  constant  k  appearing  in  (46) .  For  reliable  results, 

BK  should  be  less  than  one. 

VS(j),  ZS(j)  and  R(j)  are  arrays  that  describe  the  geometry  of  the  sur¬ 
face  S  of  revolution.  The  J  indicates  evaluation  at  t  =  >  '*cre  t  is 

the^arc  length  along  the  generating  curve,  zero  at  the  lower  pole  of  the  body 
of  revolution.  J  runs  from  1  to  2*N. 

VS  =  the  angle  (radians)  between  the  unit  vector  in  the  t  direction 
and  the  axis  {  z  axis)  of  the  body  of  revolution.  VS  is  positive 
when  is  diverging  from  the  z  axis  and  negative  otherwise. 

ZS  =  axial  distance  z  from  the  lower  pole 
R  =  distance  p  from  the  z  axis. 

All  of  the  above  data  are  printed  after  they  are  read  in.  Statement  58 

58  WRITE  (5,6U)(G(K),  K  =  1,  NG) 

prints  some  G^s  of  equation  (52).  The  Z  matrix  is  printed  by 

95  WRITE  (5,88)(Z(K),  K  =  K1,K2) 

88  FORMAT  (IX,  12Gll.lt) 

inside  DO  loops  89  and  90*  The  impedance  matrix  is  stored  columnwise  in  the 

linear  array  Z.  The  submatrix  (Z^)^  of  (45 )  is  labeled  Z1  and  printed 

columnwise.  The  first  element  of  each  column  begins  a  new  line.  For  a 

column  of  the  submatrix  (Z^).  .  the  "source  point"  j  is  fixed  and  the  "field 

n  ij 

point”  i  runs  from  the  lower  pole  to  the  upper  pole.  Similarly,  Z^,  Z^  and 


2^  are  labeled  22 ,  23,  and  Zh  and  printed.  The  impedance  matrix  is  put  into 
direct  access  storage  by 

33  WRITE  (4)(2(I),  1=1  NZ) 

Next,  the  admittance  matrix  (inverse  of  the  impedance  matrix)  is  calculated 
and  put  into  the  same  array  Z.  The  admittance  matrix  is  outputed  exactly 
as  the  impedance  matrix  was.  The  printing  is  by 

9o  WRITE  (3, 88)(Z(K),  K  =  KL,  K2) 

in  DO  loops  93  and  92.  The  direct  access  data  set  U  is  written  by 

3^  WRITE  (4)(Z(l),  I  =  1,  NZ) 

Throughout  the  program,  the  sequence 

CALL  CLOCK  (11,12) 

WRITE  (3,99)  11,12 
99  FORMAT  (IX, 217) 

appears  many  times.  CLOCK  is  a  subroutine  which,  perhaps  not  available  on 
all  360  computers,  can  be  called  as  easily  as  SIN  or  COS.  The  first  argument 
II  becomes  the  time  in  seconds  while  the  second  one  12  gives  the  year  and  day. 

A-2  Explanatory  notes  on  the  computer  program. 

LINEQ(LL,C),  the  first  subroutine  to  be  compiled,  accepts  an  LL  by  LL 
complex  matrix  stored  columnwise  in  the  linear  complex  array  C  and  returns 
the  inverse  of  this  LL  by  LL  matrix  into  the  same  array  C.  Of  course,  the  orig¬ 
inal  elements  of  the  matrix  are  destroyed.  Although  the  dunmy  array  C  is  only 


sfl 

«<  ~ 

dimensioned  C(l)  in  the  subprogram,  all  the  space  of  the  corresponding  array 
in  the  calling  sequence  in  the  main  program  is  available.  Because  of  the 
statement 

DIMENSION  LR(58) 

early  in  the  subroutine,  LL  must  not  exceed  58.  LINEQ(LL,C)  uses  the  sub¬ 
program  CABS.  Also,  since  the  matrix  elements  are  complex,  complex  addition, 
subtraction,  multiplication  and  division  routines  are  tacitly  assumed. 

The  second  subroutine  POT(AC,GS)  computes  f(0ffi)  of  (55)*  The  variables 
in  common  are 


BR1  =  d  sin  v  -  (z  -z  )  cos  v 

4.  q  '  p  q 7  q 

RR2  =-p  sin  v 
P  <1 

*  4  * <vv2 

BK  =  k 

The  arguments  of  ?0T  are 

AC  =  cos  0 
rm 

GS  =  f(0ffl) 

GS  is  the  only  output  of  the  subroutine  POT(AC,GS).  RI;1,  BR2  and  AC  are  com¬ 
bined  to  give  tQ  of  (50).  RR3  and  RR4  are  used  to  find  R  .  When  in  (53) 
tj+d  is  much  larger  than  l/1*,  the  argument  of  the  logarithm  becomes  close 
to  one.  Then,  the  series  expansion  for  log(l+x)  could  possibly  have  been 
used  in  the  subroutine  F0T(AC,GS),  but  was  not.  The  subroutine  P0l(AC,GS)  uses 
the  subprograms  ALOG,  SIN,  COS,  CMPUC,  ABS,  and  SORT. 
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First  of  all,  the  main  program  reads  the  punched  card  data  and  then 

prints  it  to  be  certain  that  it  was  entered  correctly.  The  constants  a  and  tj 

are  inserted  just  after  statement  U7.  The  DO  loop  2  puts  0.  appearing  in 

0 

(46)  into  ANG(J)  and  cos  0^  in  AC(j).  Tne  DO  loop  10  puts  ^  cos  (n-l)0^, 

j  running  from  1  to  NFHI,  in  CSM(a)  through  CSM(NIHI),  -  cor  n  0  in  CSM(NPHI+l) 

j 

through  CSM(2*NFHI),  and  ^  cos  (n+l)0.  in  CSM(2*NFHI+l)  through  CSM( 3*NPHI). 

The  nested  DO  loops  16  and  17  compute  the  Gn  of  (52).  Actually,  for  the 
fixed  n  appearing  in  (45),  Gn-1,  G  and  Gn+1  are  needed.  All  are  put  in  the 
linear  complex  array  G.  The  Gn  ^  are  put  in  G(l)  through  G(NG),  Gr  in  G(NG+l) 
through  G(2*NG),  and  Gr+1  in  G(2*£5+l)  through  G  3*NG).  Gn  is  a  KG  by 
KG  array.  (KG<=2*N  and  NG=KG*KG) .  The  complex  number  Gn  L(( J-1)*KG+I)  is 
for  the  field  point  at  I  and  the  source  in  the  vicinity  of  J.  In  DO  loop  5, 
the  subroutine  POT  is  called  upon  to  put  f(0R)  of  (53)  into  GS(K).  If  (l=J), 
the  logarithm  t cpression  of  (58)  instead  of  f(0,)  is  put  into  GS(l).  The  DO 
loop  13  performs  the  summation  (52).  The  outer  DO  loop  68  obtains  the  cases 
(n-l),  n,and  (nJl). 

After  the  G’s  are  printed  by  statement  58,  the  numbers  of  (4l)  are  in¬ 
serted  into  the  arrays  T  and  TP.  The  impedance  matrix  (45)  is  computed  in 
the  nested  DO  loops  30  and  31.  I  and  J  correspond  to  the  i  and  j  of  (42). 

The  subscripts  for  Zlt,  Z^ ,  2^, and  are  respectively  LI,  L2,  L3,  L4  so 
that  the  impedance  matrix  with  the  submatrix  arrangement  of  (25)  is  put 
columnwise  into  the  linear  complex  array  Z.  The  nested  DO  loops  70  and  71 
perform  the  double  sum  in  (4p).  II  corresponds  to  p  and  JJ  to  q. 

The  logic  between  statements  84  and  89  is  a  routine  to  print  Z.  The 


statement 


CALL  LINEQ  (NM2,Z) 


inverts  the  impedance  matrix.  The  admittance  matrix  now  occupying  Z  is  printed 
in  DO  loop  93*  The  data  included  at  the  end  of  Section  A- 3  is  for  a  sphere 
with 

NN  =  1 
N  =  10 
NFHI  =  20 

BK  =  0.39^78^2 

BK  is  the  propagation  constant  for  a  radius  of  0.2  wavelengths.  The  portion 
of  the  data  defining  the  geometry  of  the  sphere  was  actually  computed  and 
punched  on  cards  by  a  short  auxiliary  computer  program. 

If  N  or  NFHI  is  too  large,  the  dimension  of  sone  arrays'' must  be  increased. 
All  the  arrays  whose  dimension  depends  upon  N  or  NFHI  are  listed  below. 

COMPLEX  Z(NM2**2),  GS(NFHI),  G(l2*N*N) 

DIMENSION  VS(N*2),  SV(N*2),  CV(N*2),  ZS(N*2), 

R(N*2),  ANG(NFHI),  AC(NFHI-), 

CSM(3*NFHI) 


Here,  as  in  the  program, 
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CALL  LINEQ  (NM2,Z) 

inverts  the  impedance  matrix.  The  admittance  matrix  now  occupying  Z  is  printed 
in  DO  loop  93-  The  data  included  at  the  end  of  Section  A- 3  is  for  a  sphere 
with 

NN  =  l 
N  =  10 
NFHI  =  20 
BK  =  0.39^78^2 

BK  is  the  propagation  constant  for  a  radius  of  0.2  wavelengths.  The  portion 
of  the  data  defining  the  geometry  of  the  sphere  was  actually  computed  and 
punched  on  cards  by  a  short  auxiliary  computer  program. 

If  N  or  NFHI  is  too  large,  the  dimension  of  sone  ari ays"  must  be  increased. 
All  the  arrays  whose  dimension  depends  upon  N  or  NPHI  are  listed  below. 

COMPLEX  Z(NM2**2),  GS(NFHI),  0(l2*N*N) 

DIMENSION  VS(N*2),  SV(N*2),  CV(N*2),  ZS(N*2), 

R(N*2),  AMG(NFHI),  ac(nfhi-), 

CSM(3*NFHI) 


Here,  as  in  the  program, 


NM2  =  2*(N-1). 


I'iK-xnaj 


[  A-3.  The  computer  program 

I  //REV170  JOB  (0034, EE,  1 1 ,03 >  ,  *MAUT? ,  JOF  •  .MSr.LFVFl*  1 

|  //  EXEC  FORTGCLO  *  P  ARM.  FOR  T= 'MAP  •  ,  P  ARM,  |_KF0=  '  X&EF  • 

t  / /FORT ,SY$In  00  * 

SURROUT  I  ME  UNFQ(LL.C) 

COMPLfcx  ell), STOP, STD, ST, S 
DIMENSION  LR(5fi) 

00  20  I  tel  ,LL 

LRU)*! 

20  CONTINUE 
M 1  *0 

00  18  M»1 ,LL 
K  =  M 

00  2  I=M,U 
K 1  ssM  1  ♦  1 
K2e=MH-K 

!E(CABS(C(K1 M-CAHS(C(X2) » >  2,2,4 

,  (S  K  =  I 

2  CONTINUE 
LS=LR(M) 

*  LR(M)=LR(K) 

LR(K)=LS 

K2=M1+K 

ST0R=C(K2) 

.1 1  *0 

on  7  jet i  ,ll 
ki=ji+k 

K2=J1+M 

ST0=C(K1» 

C ( K 1 ) *C ( K? ) 

C(K?)=STO/STU« 

J1*J1+LL 
7  CONTINUE 
K1=M1+M 
C(K1»=1,/STUR 

*  CO  11  1*1  ,  L  L 

IE(I-M)  12,11,12 
12  M*M1  +  I 

!•  ST*C(K1» 

C(K1 )*0. 

J1=0 

00  10  J*1  ,LL 

KlaJl+I 

K2rJl+|# 

C(Kl)*C(Kl)-C(K2)*ST 

J1SJ1*LL 

10  CONTINUE 

11  CONTINUF 
Ml sMI+LL 

18  CONTINUE 
J 1  *0 

>  on  9  J=1,LL 

IE ( J— LR ( J I )  14,8, j  & 

14  IRJ«IR(J) 

A  J2*( LR J-l  )*LL 

21  no  13  l»l,Ll 
K2*J2*l 
K1*J1+I 
S«C(K2) 

C I K2 ) *C ( K 1 ) 

C(Kl  1*S 


t 


>•<»"' -’r'tr'Kt 
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13  CONTINUE 

LR<  J)«LR(LRJ) 

LR(LRJ)*LRJ 
IF  ( J-LR IJI!  14,B,14 
8  Jl® Jl+LL 
R  CONTINUE 
RETURN 
END 

SUBROUTINE  POTUC,GS) 

COMMON  RR1,RR2,RR3,RR4.?K 
COMPLEX  W3,W4,6S,G 
YS«RR1+AC*RR2  \ 

Y=A8S ( YS I 
R0*RR3+RR4*AC 
RK*BK*SORT (RD) 

02*RO-Y*Y 
Yl*Y-.25 
Y2»Y+.25 

Rl*SCRT(Yl*Yi+02) 

R2*SORT(Y2*Y2+D2) 

IF(Yl)  6,7,7 

A  TIN*2.*ALOG( (-Y1+R1 ) *< Y2+R2 I /D2 ) 

GO  TO  R 

7  TIN  ®2.*ALOG< < Y2+R2  ) / ( Y 1+R 1 ) ) 
t  8  SN«SIN(RK) 

!  CS*COS( RK ) 

)  S*BK-RK*T IN 

\  W3«CMPLX<CS,-SN) 

I  W4«CMPLX(TIN,-S) 

{  GS»W3*W4 

j  RETURN 

I  END 

i  COMMON  RR1,RR2,RR3*RR4,BK 

COMPLEX  A3,A4,Z!3364),GS(40),G(4«00) 

DIMENS IONVS! 60 ),SV! 60)  ,CV(  AO)  ,ZSC  AO)  ,R  1 60)  ,  ANG  ( 40 )  *ACUO)  ,CSM(  120) 
I  DIMENSION  TPI4) , T ( 4 ) , JK (4 ) 

1  l  REWIND  4 

|  50  READ! 1 , 51  * ENO=52 )NN,N, NPHI , BK 

?  51  FORMAT (3I3»E14,7) 

I  KG»2*N 

61  READ(1,53)(VS(J).J*1,KG) 

|  62  R£A0(l,53)fZS(J),J«l,KG) 

!  63  READ! 1,53)(R(J),J*1,KG) 

53  FORMAT ( 10F8.4) 

WRITE(3,54)NN,N,NPHI,BK 

54  FORMAT ( IX// •  NN«',I3,‘  N«*,I3.*  NPHI«*,I3»*  BK**,E14.7) 

55  FORMAT ( IX/ *  VS‘) 

56  FORMAT  I IX/*  ZS») 

|  57  FORMAT! IX/*  R*) 

.  WRITE(3,55) 

WRITE (3,46) ( VS ( J ) , J*1 ,KG ) 

46  FORMAT (1X,10FB«4) 

WRITEI3,56) 

WRITE(3,46)(ZS(J),J-1,KG) 

WRITE(3»57) 

WR|TEI3,46)(R(J),J«1,KG) 

DO  47  J*1,KG 
SVi J )«SIN( VS ( J ) ) 
i  CV(J)«COS(VS(J)) 

47  CONTINUE 


PI»3. 141593 

ETA-376. 707 

NM-N-1 

NM2*NM*2 

NZ«NM2*NM2 

NG«KG*KG 

M5-NN+2 

M6-NN+4 

FM-NN 

FM2«NN*NN 

DP-PI/NPHI 

00  2  J»1,NPHI 

ANG( J ) * ( J-.5 ) *0P 

AC! J)«C0S!ANG! J) ) 

2  CONTINUE 
M3-0 

00  10  MM-M5,M6 
Ml-HM-3 
M2«M3*NPH1 
DO  11  K* 1 tNPHI 
K1-M2+K 

CSMIKi)-OP*COS<Ml*ANG!K)) 

11  CONTINUF 
M3-M3+1 
10  CONTINUE 

CALI  CLOCK! 11*12) 

WRITE!3, 99)11. 12 
99  FORMAT ( IX, 217) 

00  16  J-l.KG 
AA-DP*R!J)*,5 
OG  17  I-l.KG 
Z3»ZS!J)-7S!I> 

RR1-5V!  JI'*R(  J)+CV!  J)*Z3 
RR2— 'V!  J)*R!  I ) 

RR3*R( J)*Rt J)+R( I )*R< I )+Z3*Z3 
RR4— ?.*R!  J)*R! I ) 

21  DO  5  K-2.NPHI 

CALL  POT!AC(K),GS(K)> 

5  CONTINUE 

IF(I-J)  7,8,7 
8  X-R! J )*DP 

XX«S0RT!l./16.+X*X) 

W1«(2.*X*AL0G( (,25+XX)/X)+.5*AL0G(!X+XX)/.25))/AA 
W2«-BK 

GS(  1  )«W1-M0.,1.)*W2 
GO  TO  67 

7  CALL  POT ! AC! 1 ) ,GS( l ) ) 

67  M3«! J-1)*KG+I 
DO  68  MM«l,3 
Ml-MM-1 
M4«M1*NPHJ 
M2«M1*NG+M3 
G!M2)«0, 

DO  13  K-l.NPHI 
K2«K+H4 

G ! M2 ) *G ( M2 ) +GS ! K ) *CSM I K  2 ) 

13  CONTINUE 

68  CONTINUE 
17  CONTINUE 
16  CONTINUE 


CALL  CLOCK!! 1* 12) 

WRITE13,  99)11, 12 
WRITE! 3 ,48) 

48  FORMAT ( IX/  •  POTENTIALS  FROM  THE  CIRCULAR  RIBBONS  OF  CHARGE*  ) 

58  WRITE<3,64) <G!K) ,K*l,NG) 

64  F0RMAT(1X/(1X,10G11.4)) 

CALL  CLOCK! 11,12) 

WRITEI3, 99)11, 12 

TP! 1 ) *.5 

TP!2)*.5 

TP<3)*-.5 

TP!4)*-,5 

T!l)*.125 

T!2)*.375 

TI3)«.375 

T ( 4)*« 125 

CA*BK*ETA 

CO-CA/BK/BK 

DO  30  J*1,NM 

J1*2*!J-1) 

JL*! J-1)*NM2 

00  31  I«1,NM 

L1*JL+I 

L2-L1+NM 

L3*L1+NM2*NM 

L4*L3+NM 


Z!L1)*0. 

2!L2)*0. 

Z!L3)*0# 

Z!L4J*0. 

I1*2*!I-1) 

DO  70  JJ-1,4 
J2-J1+JJ 
KA1«|J2-1)*KG 
KB1*KA1+NG 
KC1*KB1+NG 
DO  71  11*1,4 
12*11+1! 

SS*SV! I2)*SV( J2) 

CC*CV( !2>*CV!J2) 

KA2*KA1+I2 
KB2-KB1+I2 
KjC2*KC1  +  I2 

A3*.5*!6IKC2)+G!KA2)) 

A4*.5«!G!KC2)~G!KA2) ) 

74  Z(L1)*Z(L1 )  +  (CA*T ! 1 1 )*T ! JJ )*( SS*A3+CC*G!KB2 ) )-CO*TP! II )*TP! JJ  )*G!K 
1B2))*!0.,1.) 

Z!L2)«Z(L2)+CA*SV!J2)*T!II )*T! JJ )*A4-FM*C0*G< KB2 >*T » 1 1 )*TP!JJ)/R! I 
12) 

Z(L3I*Z!L3)-CA*SV!I2)*T!I1 )*T!JJ)*A4+FM*C0*G!KB2)*TP( I! )+T! JJ)/RI J 
121 

Z!L4)«Z!L4)+1CA*A3-FM2*CQ/R! I2)/R< J2)*G1KB2) )*T! !I)*T! JJ)*iO., 1. ) 
71  CONTINUE 
70  CONTINUE 
31  CONTINUE 
30  CONTINUE 

CALL  CLOCK! 11,12) 

WRITEI3, 99)11, 12 
84  JK! 1 )*1 
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JK!3)*NM2*NM+1 
JK!4)*JK(3)+NM 
on  89  J*1.4 
K1*JK< J) 

WRITE<3,49) J 

49  FORMAT !  IX/ •  Z‘,I1) 

DO  90  1*1, NM 
K2*K1+NM-1 

95  WRITE(3»P8)(Z!K),K*K1.K2) 

88  FORMAT (1X« 126 11*4) 

K1»K1+NM2 

90  CONTINUE 

89  CONTINUE 

33  WRITE  !4)!Z!I),I*1,NZ) 

CALL  CLOCK ( 1 1  * !2 ) 

WR*Te(3,99)  11,12 
CALL  LINE0(NM2,Z) 

CALL  CLOCK! 11,12) 

WR I TE ( 3 ,99 )  11,12 
DO  93  J*1 ,4 
K1«JK( J) 

WRITE! 3,24) J 

24  FORMAT ! IX/ *  Y*,I1) 

00  92  1*1 ,NM 
K2*K1+NM-1 

96  WRITE!3,88)!Z!K),K*K1,K2) 

K1-K1+NM2 

92  CONTINUE 

93  CONTINUE 

34  WRITE  !4MZm,I*l,NZ) 

CALL  CLOCK! 11,12) 

WRITEI3, 99)11, 12 

GO  TO  50 

52  STOP 
END 

/* 

//G0.FT04P001  00  DSNAME*EE0034.REV,DISP*0LD,UNIT*2314,  X 

//  V0LUME*SER*SU0004,DCB*!RECFM«V,BLKSIZE*1800,LRECL*1796) 

//GO.SYSIN  00  * 

001010020  0.3947842E+00 


1.4923 

1.3352 

1.1781 

1.0210 

0.8639 

0.7069 

0.5498 

0.3927 

0.2356 
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0.0785 
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4.4012 
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5.2504 

5.6035 

5.8971 
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0.2497 
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1.2181 
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2.4204 
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v  The  problem  of  electromagnetic  radiation  and  scattering  from  perfectly  con¬ 
ducting  bodies  of  revolution  of  arbitrary  shape  is  considered.  The  mathematical 
formulation  is  an  integro-differential  equation,  obtained  from  the  potential 
integrals  plus  boundary  conditions  at  the  body.  A  solution  is  effected  by  the 
method  of  moments,  and  the  results  are  expressed  in  terms  of  generalized  networjc 
parameters.  A  computer  program  for  computing  the  generalized  impedance  matrix  of 
an  arbitrary  body  of  revolution  is  included. 

The  expansion  functions  chosen  for  a'  moment  solutions  are  harmonic  in  0 
(azimuth  angle)  and  subsectional  in  t  (contour  length  variable ) .  Beeause  of 
rotational  symmetry,  the  solution  becomes  a  Fourier  series  in  0,  each  term  of 
which  is  uncoupled  to  every  other  term.  Hence,  the  problem  reduces  to  a  set 
of  independent  modes,  one  for  each  harmonic  term. 

Illustrative  computations  are  given  for  radiation  from  apertures  and  plane- 
wave  scattering  from  bodies  of  revolution.  The  impedance  elements,  currents, 
radiation  patterns,  and  scattering  patterns  for  a  conducting  sphere  are  computed 
both  front  the  general  program  and  from  the  classical  eigenfunction  solution.  The 
agreement  obtained  serves  to  check  the  general  program.  Similar  computations  for 
a  cone-sphere  illustrate  the  application  of  the  general  program  to  problems  not 
solvable  by  classical  methods. 
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